Shape of Heteroepitaxial Island Determined by Asymmetric Detachment 
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Square lattice gas models for heteroepitaxial growth are studied by means of kinetic Monte 
Carlo simulations, in order to find a possible origin of anisotropic island shape observed in growth 
experiments of long organic molecules. When deposited molecules form clusters irreversibly at 
their encounter during surface diffusion, islands grow in a ramified dendritic shape, similar to 
DLA. Introduction of molecular detachment from edges makes islands compact with smooth edges. 
Tilting of adsorbed long molecules or steps in a vicinal substrate may induce orientation-dependence 
in the detachment rate of edge molecules from an island. In simulations with orientation-dependent 
detachment rates, a clear anisotropy in an island shape is observed. Shape anisotropy on a vicinal 
substrate is enhanced as steps get dense, in agreement to the experimental observation. 



PACS numbers: 81.15.Aa, 68.43.Hn, 68.47.Pe 



I. INTRODUCTION 

There have been increasing interests in heteroepitaxial 
growth of organic semiconductors in relation to the fab- 
rication of thin and deformable microelectronic or opto- 
electronic devices [l|,|2|,|3j. To accomplish a high standard 
in electrical properties, improvement of the film quality 
is indispensable, and fundamental researches on growth 
mechanisms of adsorbed organic molecules are under- 
taken [4j,l5[. In an experiment of pentacene (Pn, C14H22) 
growth on Si(OOl) [4(, authors concluded that an organic 
thin-film growth is similar to the epitaxial growth of in- 
organic materials, such as the formation of fractal rami- 
fied islands like diffusion-limited aggregates (DLA) @, Q . 
On the other hand, in an experiment of Pn growth on a 
hydrogen-terminated Si(l 11) surface [H-Si(lll)] [5[, ape- 
culiarity is found. On a fiat H-Si(lll) surface islands are 
compact with smooth edges and isotropic, whereas on 
a vicinal surface islands are dendritic in the step-down 
orientation while edges in the step-up orientation remain 
compact: Each island has a strong anisotropy in its shape 
on a vicinal substrate. Island shape anisotropy, however, 
is not specific to a vicinal surface. It is observed in the 
vapor growth of another organic molecule, a behenic acid, 
even on a fiat singular substrate, namely on an oxidized 
GaAs Q. In the present paper we analyze some model 
systems which lead to shape anisotropy theoretically. 

In the heteroepitaxial growth of inorganic materials, 
island shape is mainly controlled by the competition of 
two processes after an adatom deposition; an adatom dif- 
fusion on a terrace surface and an incorporation kinetics 
at island edges 0, [TO, EH, EH • With a fast incorporation 
at island edges, diffusion controls the growth and edges 
undergo a morphological instability to form dendrites. In 
the extreme case of irreversible solidification such that an 
adatom once solidified never detaches again from the is- 
land edge, then an island takes a ramified irregular shape 
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with many branches pjjj]. Its structure is similar to the 
DLA without a characteristic length, and is called a frac- 
tal. An island is isotropic or symmetric by reflecting the 
lattice symmetry. 



To provide a compact island shape, edge smoothening 
processes are necessary, such as an edge diffusion or a de- 
tachment from the edge. If the detachment rate is high, 
adatoms loosely attached to the edge will readily de- 
tach and they are hardly incorporated in two-dimensional 
(2D) islands. Then, the slow incorporation kinetics gov- 
erns the growth to make islands compact with smooth 
edges [li] Il5l [l6| . Islands may be round at high temper- 
atures or polygonal at low temperatures. 



In the heteroepitaxial growth of organic molecules, one 
has to take an additional feature into account, i.e. the 
size of an adsorbed organic molecule. It can no longer 
be regarded as a point object as in inorganic cases so far 
discussed. Both pentacene and behenic acid molecules 
are fiat and elongated in shape. When they are adsorbed 
on a substrate, molecules can be normal or lateral to 
the substrate surface [18[. Pn molecules are known to 
grow a wetting film with the molecular long axis nor- 
mal to the H-Si(lll) surface, namely in a standing- up 
orientation [l7[. Behenic acid is also known to be in a 
stand-up orientation but with a small tilting on a sub- 
strate. We assume that this molecular orientation affects 
the detachment rate of molecules from island edges, and 
leads to anisotropy in an island shape. 

In a previous study [l^ we proposed a simple lattice 
gas model, and studied general aspects of the effect of 
anisotropic detachment on an island shape. In the next 
§2 we introduce the model with simulation results of is- 
land shapes on a fiat substrate. The model is extended 
in §3 so as to include the effect of steps explicitly for the 
case of adsorption on a vicinal substrate. The last §4 
summarizes the result. 
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II. 2D ISLANDS ON A FLAT SUBSTRATE 

When molecules in a 2D crystal stack almost normal 
to a substrate surface but with a finite tilting, as to the 
a-axis in the case of a behenic acid [8j], kinetics at island 
edges may depend on relative orientations of the edge and 
of the molecular tilting. To an edge where the molecu- 
lar inclination restricts the incorporation space, diffusing 
molecules may be difficult to be attached to the edge, 
but once attached they may be difficult to be detached 
from this edge. Depression of the attachment rate slows 
down the incorporation kinetics, and the edge becomes 
smooth. On the contrary, suppression of detachment in- 
creases incorporation and leads to a diffusional instability 
of a smooth edge. So far the relation between the molec- 
ular tilting and the orientation of dendritic edge is not 
identified experimentally [8[ , we here assume simply that 
a molecular tilting affects detachment process such that 
the detachment rate of a molecule depends on the edge 
orientation. In our theoretical treatment, we forget about 
the finite size and shape of an adsorbing organic molecule, 
but consider it implicitly in the orientation-dependence 
of the detachment rate. 

We perform kinetic Monte Carlo simulations pjl [2(j 
of a lattice gas models of point molecules depositing on a 
square substrate with a deposition rate F per area. Ad- 
sorbed molecules diffuse on a surface with a diffusion con- 
stant D s . When two diffusing molecules meet, they make 
a bond to lower an energy and form a cluster. Clusters 
enlarge their sizes by incorporating diffusing molecules 
at edges. On the other hand, molecules at edges can 
detach from the islands when they are loosely attached: 
Those particles with a single bond to the island can break 
the bond to migrate out on the substrate surface with a 
rate D ei but molecules with more than two bonds are 
assumed immobile. When a particle is deposited above 
an island, it diffuses on the island terrace till it reaches 
an edge, and steps it down to be incorporated into the 
island. Evaporation of adsorbed molecules from the sur- 
face is excluded, since the growth is assumed to take 
place at a low temperature. When there is no detach- 
ment, grown islands take a ramified dendritic form with 
irregular fine branches, similar to the DLA [6|,l3|. With a 
small detachment, one obtains dendritic aggregates with 
fat branches, as shown in Fig. 1(a). Here, parameters 
are D s /F = 10 12 and D e /F = 10 4 with a system size 
L 2 = 500 2 and the coverage O = 0.1. Although islands 
interfere with each other in a diffusion field, they have al- 
most a symmetric shape. By applying the box-counting 
method, one obtains a rough estimate of the fractal di- 
mension as df « 1.70. 

We now introduce an anisotropy in the rate of detach- 
ment from edges. On the left edge of an island, for in- 
stance, we arbitrary reduce the rate of detachment by Rl- 
Island morphology simulated by means of kinetic Monte 
Carlo simulation [20[ becomes anisotropic, as shown in 
Fig. 1(b) and 1(c). When the detachment from the 
left edge is reduced by a factor 10 as Rl = 0.1 (Fig. 
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FIG. 1: Islands on a substrate with a size 500 2 with a surface 
diffusion D s /F — 10 12 at a coverage = 0.1. A detach- 
ment rate is D e /F = 10 4 , and (a) isotropic, and (b) and (c) 
anisotropic. Detachment from the left edge is suppressed by 
a factor (b) R L = 0.1, and (c) R L = 0.01. 



1(b)), the left edges of islands grow faster than those in 
the other orientations, and the primary branch to the 
left extends about three times more than to the right. 
The fractal dimension by box counting remains the same 
value df « 1.7. On further reducing the detachment 
rate from the left edges to Rl — 0.01 (Fig. 1(c)), is- 
lands point sharply to the left with the left-to-right ratio 
ranging about 3 to 6, depending on environment. Pri- 
mary branches extending to the left are fine and have 
less side-branches, compared to Fig. 1(a) and 1(b), and 
primary branches to other orientations have many sec- 
ondary branches extending long to the left. Pointing 
structure leads to the small fractal dimension of df « 1.6 
by the box-counting method. Island morphology in Fig. 
1(c) looks quite similar to those observed for the behenic 
acid islands on an oxidized GaAs [8[. 

Instead of suppressing detachment from the left edge, 
the enhanced detachment from the righ edge leads to 
anisotropic island shape [19[. In this case, however, is- 
lands extend vertically in y direction, and the obtained 
shape looks different to what is observed in the exper- 
iment [8]. On increasing the detachment rate D e /F, 
islands become more compact with smooth edges as a 
skeletal shape or even to a complete square. The effect 
of anisotropic detachment in these cases have been dis- 
cussed previously [l^ . 



III. ISLANDS ON A VICINAL SURFACE 

We know consider the case of a heteroepitaxial growth 
on a vicinal surface, which might be relevant to Pn 
growth on H-Si(lll) 0|. On a flat terrace, Pn molecules 
form 2D islands, and they are compact with an isotropic 
shape. Therefore, Pn molecules might often detach from 
island edges. On a vicinal H-Si(lll) surface, Pn islands 
show different morphology in the upward and downward 
directions relative to the substrate steps. Island edges 
extending to the step-down direction are pointed in a 
dendritic form, while those to the step-up directions re- 
main to be smooth and round. 

Since Pn is known to grow in a standing-up orientation 
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FIG. 2: Detachment rates of adsorbed long molecule from 
islands on a vicinal surface. 



[17l | , an interaction between Pn molecules is stronger than 
that between Pn and the substrate. Let us now imagine 
that a Pn crystalline film is growing on a vicinal surface 
with a positive slope, as shown in Fig. 2. When a film 
grows to the left in a step-down direction and reaches an 
upper side of the descending step, as a left island in Fig. 
2, further crystal growth takes place by incorporating a 
Pn molecule from the lower terrace. In this case, the 
whole length of a newly attached molecule interacts with 
the island edge and the step ledge, and it is hard to be 
detached back to the lower terrace. Also, the free upper 
part of the long molecule may in general be susceptible 
to thermal fluctuation, and there may be a substantial 
entropy contribution which reduces the bonding free en- 
ergy. Attachment from the lower terrace may lessen this 
entropy contribution and enhances the Pn-Pn bonding. 
Thus, the detachment rate from the lower terrace is di- 
minished from the normal rate D e on a flat terrace by 
a factor Rl < 1 as to D g Rl. On the contrary, when 
the island grows to the right in the step-up direction and 
reaches the lower side of the ascending step, as a right 
island in Fig. 2, the next Pn molecule to be incorporated 
is located on the higher terrace, and its upper part does 
not contribute to the Pn-Pn bonding. Since the substrate 
step has a height of about one fourth of the length of a Pn 
molecule, a loss in the bonding energy may be significant. 
Accordingly, those molecules incorporated into an island 
edge from the higher terrace are easily detached back to 
the higher terrace with an enhancement factor Rh > 1. 
We now study the effect of this anisotropic detachment 
with rate modification factors Rh and Rl on the island 
morphology. 

In kinetic Monte Carlo simulation of heteroepitaxial 
growth on a vicinal surface, steps are separated with a 
distance L s . After a deposition, molecules perform D s /F 
times of diffusive migrations until a monolayer is covered, 
but during this period molecules collide to form clusters. 
Steps are assumed to have no effect on the surface dif- 
fusion, for simplicity. A large value of surface diffusion 
constant D s /F = 10 15 is chosen in order to grow only 
a single island in our system of a size L 2 = 500 2 . From 
edges of an island on a flat terrace, molecules with a single 
bond detach with a rate D e /F — 10 5 . Above a descend- 
ing step, detachment is enhanced by a factor Rh = 10, 
but in front of the step detachment is assumed completely 
forbidden, Rl = 0, corresponding to an extreme case. 

By initially providing an embryo of a size 2 by 2 at the 
center of the system, an island is nucleated and grows to 
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FIG. 3: An island on a vicinal substrate surface. The surface 
diffusion is D s /F — 10 15 and the detachment rate on a flat 
terrace is D e /F — 10 5 . Detachment to the higher terrace (to 
the right) is enhanced by a factor Rh — 10 and to the lower 
terrace (to the left) is completely suppressed: Rl — 0. Step 
separations are L s =(a) 5, (b) 4, (c) 3 and (d) 2. The coverage 
is B = 0.1 on a substrate with a size 500 2 . 



the coverage 6 = 0.1. As the step separation decreases 
from L s = 5 to 2, the island shape alters as shown in Fig. 
3(a) to 3(d). Even with a very large surface diffusion, a 
second island is nucleated accidentally in Fig. 3(c). At a 
large step separations L s = 5 and 4, the island has only 
a weak left-right anisotropy, as in Fig. 3(a) and 3(b). 
When steps come closer as L s = 3 and 2, anisotropy 
becomes obvious, as in Fig. 3(c) and 3(d). 




FIG. 4: Ratio of the edge velocities to the lower side Vl to 
the higher side Vh as a function of the step separation L s . A 
curve is a guide for eyes. The surface diffusion is D s /F — 10 15 
and the detachment rate is D e /F — 10 5 with an enhancement 
factor to the uphill Rh = 10 and a complete suppression to 
the downhill Rl — 0. 

In order to quantify the shape anisotropy, we measure 
the ratio Vl/Vh of the front velocities to the lower side 
Vl and to the higher side Vh as a function of the step sep- 
aration. The result shown in Fig. 4 clearly demonstrates 
that an anisotropy increases as the step separation de- 
creases, in a qualitative agreement to the experimental 
result obtained in Pn system [5[. 

So far, we discussed the shape of a single island. In 
a large area of the vicinal surface many islands are nu- 
cleated and they affect their mutual growth and shape 
through the diffusion field. In the Pn experiment one 
may identify some order in spatial arrangement of islands 
[5[ : Islands look to align in queues parallel to the tilting 
direction of a vicinal substrate, with a dendritic left edge 
of one island being contiguous to the smooth right edge 
of a left one. ( See, for example, Fig. 2 (c) and (d) in 
Ref.0. ) 
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(a) (b) (c) 

FIG. 5: Island distribution on a vicinal surface with a size 
500 2 at a coverage G = 0.3. Step separation is L s — 2 and 
surface diffusion constants are (a) D s /F = 10 10 , (b) 10 11 , and 
(c) 10 12 . Other parameters are D e /F = 10 4 , R H = b,R L = 
0. 

In our system with a limited size L 2 = 500 2 , island den- 
sity is controlled by varying the surface diffusion constant 
D s . In some combinations of parameter values, simula- 
tion can produce island arrangements with some spatial 
ordering, as shown in Fig. 5. The height of the vicinal 
surface is increasing to the right with ascending steps at 
the separation L s = 2. Molecules are deposited to the 
coverage 6 = 0.3. The surface diffusion constant is set at 
(a) D s /F = 10 10 , (b) 10 11 , and (c) 10 12 , by keeping the 
detachment rate fixed to D e /F — 10 4 , and the enhance- 
ment factor Rh — 5 and no detachment to the step-down 
direction Rl = 0. At the smallest surface diffusion D s /F 
in Fig. 5(a), there are many islands distributed almost 
randomly. The stochastic noise associated to random de- 
position is frozen in during the nucleation process. In 
Fig. 5(b), as the surface diffusion increases, island den- 
sity decreases and the diffusion process seems to average 
out the deposition noise and to induce some correlation 
in islands arrangement. Islands do not seem to align per- 
pendicular to the step orientation. One might notify sim- 
ilar tilting in islands arrangement in the Pn experiment 
0| . For larger D s in Fig. 5(c), island density becomes so 
low that it is difficult to identify spatial correlation for 
certain. 
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(a) (b) (C) 

FIG. 6: Spatial correlation function g(r) corresponding to 
Fig. 5. "H" marks the high hill-top, and "L" marks the low 
valley. 

In order to evaluate spatial correlation explicitly, we 
calculated the correlation function 

g(r) = (£n(r + r i )n(r i ))/N (1) 



where the occupation variable n(r^) is unity when the lat- 
tice site i is occupied by an admolecule, and zero when 
the site i is empty. Since the system is periodic both 
in x and y directions with a periodicity L, the correla- 
tion function has the symmetry g{x,y) = g(L — x,y) = 
g(x,L—y) = g(L—x,L—y). Thus, we obtained the height 
contour of g(x,y) for < x,y < L/2, as shown in Fig. 
6. A mark "H" represents the high hill-top with a strong 
correlation, and a mark "L" represents the low valley 
with a weak correlation. From Fig. 6(b) one clearly ob- 
serves that the pattern has a strong periodicity with a pe- 
riod p = (L/4,L/8) at a surface diffusion D s /F = 10 11 . 
With a smaller diffusion D s /F = 10 10 , average island 
density gets higher as in Fig. 5(a), and the periodicity 
is not commensurate to the system size. Therefore, g(r) 
has a rather broad first peak around r = (100, 100) in 
Fig. 6(a). With a larger diffusion D s /F = 10 12 , there 
are a few islands as in Fig. 5(c), and g(r) has a broad 
peak at around r = (150, 200). 

So far, we discussed an extreme case with the step sep- 
aration L s = 2 in order to realize a strong anisotropic 
effect in our small system. Since the island shape 
anisotropy is expected to be caused by an accumula- 
tion of many weak anisotropy effect at steps, one has 
to pack as many steps in an island as possible. Actually, 
with L s = 3 we can observe similar arrangement of many 
anisotropic islands, but with a larger separation L s > 4, 
it becomes harder to observe anisotropy in islands. In a 
real system of Pn molecules, step separations are wider, 
4.5nm to lOnm, but an island size is about 10/im and 
it spans many steps. It may be interesting if there is a 
critical size for an island to show a shape anisotropy. 



IV. SUMMARY AND DISCUSSIONS 

Motivated by heteroepitaxial growth experiments of 
long organic molecules on a flat [8] or on a vicinal surface 
[5[, effect of detachment anisotropy on the island shape 
is studied by means of kinetic Monte Carlo simulations 
of a lattice gas model. 

Only with a molecular deposition followed by a surface 
diffusion, islands take a ramified dendritic shape, similar 
to DLA. Addition of detachment from the edge makes fat 
dendritic branches, and with a large enough detachment 
rate islands become compact in a square form. 

If the long molecules are tilted in an epitaxial film, the 
rate of detachment from island edges may depend on the 
edge orientation. Simulation shows that the edge with a 
smaller rate of detachment shows diffusional instability 
and forms a dendritic tip. Resulting island shape looks 
similar to the one obtained in the experiment [8]. 

On a vicinal surface, an interaction between an edge 
molecule and a neighboring one depends on edge position. 
When an edge has just climbed up the substrate step, a 
molecular bonding is weak and an edge molecule may 
easily detach from the edge site. On the other hand, 
when an edge has just stepped down the substrate step, 
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a bonding is strong and the edge molecule is hard to 
detach. In our simulation we assume an extreme case of 
no detachment at a lower edge to stress an anisotropy 
effect. Simulations show clear anisotropy in adsorbed 
island shape as the step density increases. 

In addition to the shape anisotropy of a single island, 
the orientation-dependent detachment rate may induce 
a spatial correlation in the arrangement of islands with 
clear peaks in the spatial correlation function. Since the 
random deposition process introduces shot noise in the 
nucleation process, a fast surface diffusion seems neces- 
sary to average out fluctuation and to induce a spatial 
order, but then the island separation becomes large. For 
further studies on spatial correlation, one needs a larger 
system size. 

As for the reason of island shape anisotropy, one may 
think about a diffusional anisotropy, such that a step 
provides an asymmetry in an energy barrier for the ad- 
molecule diffusion. However, the asymmetric potential 
is unable to lead unidirectional motion nor a net drift, 
as the Feynman's ratchet and pawl machine cannot pro- 



voke directional rotation [2l|. Another effect that steps 
may provide is a nucleation center by lowering the energy 
barrier. But in the experiment there is no observation of 
alignment of nucleation centers along the step direction. 
Therefore, the latter effect might be small. Observation 
of admolecule motion on the substrate, if possible, should 
resolve these points clearly. 
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